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Abstract 

We have performed first-principles calculations of graphene edge stresses, which display two 
interesting quantum manifestations absent from the classical interpretation: the armchair edge 
stress oscillates with nanoribbon width and the zigzag edge stress is noticeably reduced by spin 
polarization. Such quantum stress effects in turn manifest in mechanical edge twisting and warping 
instability, showing features not to be captured by empirical potentials or continuum theory. Edge 
adsorption of H and Stone- Wales reconstruction are shown to provide alternative mechanisms in 
relieving the edge compression and hence to stabilize the planar edge structure. 
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Graphene, a two-dimensional (2D) single layer of carbon atoms, has attracted tremen- 
dous attention because of its nni q ne electron, propertiesfl and potential applications in 
electronic devices [2J. Earlier studies have focused on characterizing the unusual electronic 
and transport properties of graphene, particularly as a massless Dirac fermion systemjl, 2]. 
More recently, some attention has been shifted to the structural stability of graphene [2j-l5|, 
which is critically important to realizing the potential applications of graphene. On the one 
hand, as a 2D membrane structure, graphene provides an ideal testing ground (3), Q for the 
classica. Mernrin-Wagner theorem on the e* 1S te„ce of .ong-range crysta.line order in vM. 
On the other hand, the free edges of graphene are amenable to edge instabilities [5|, I8H10|. 

The graphene edge stability is characterized by two fundamental thermodynamic quan- 
tities: edge energy and edge stress. The edge of a 2D structure can be understood in 



analogy to the surface of a 3D structure 



11 



12]: the edge (surface) energy accounting for 



the energy cost to create an edge (surface) defines the edge (surface) chemical stability; 
the edge (surface) stress accounting for the energy cost to deform an edge (surface) defines 
the edge (surface) mechanical stability. First-principles calculations showed that chemically 
the armchair edge is more stable with a lower energy, while the zigzag edge is metasta 
against reconstruction Q], and both edges are hydrogenated in an H-rich environment 
Empirical-potential calculations showed that both intrinsic edges are under compressive 



io(. 



stress rendering a mechanical edge twisting and warping instability 

Usually, stress and mechanical instability are understood as phenomena of classical me- 
chanics, but they are expected to be affected by quantum effects which become prominent 
at nanoscale. So far, however, quantum effects have been mostly shown for electronic struc- 
ture and energetic quantities of low-dimensional nanostructures. Here, we demonstrate an 
interesting example of quantum manifestations of mechanical quantities in graphene edge 
stress. Using first-principles calculations, we predict that the armchair edge stress in a 
nanoribbon exhibits a large oscillation with ribbon width arising from quantum size effect, 
while the zigzag edge stress is reduced by spin polarization. Such quantum effects on edge 
stress in turn manifest in graphene edge mechanical instability, with "quantum" features 
that apparently cannot be described by empirical force-field potentials or continuum theory. 

Our calculations were performed using the method based on density functional theory 
(DFT) in the generalized gradient approximation, with the Perdew-Burke-Ernzerhof func- 
tional for electron exchange and correlation potetnials, as implemented in the VASP code 
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FIG. 1: The armchair edge stresses and edge energies of graphene nanoribbons as a function of 
ribbon width. Inset: schematics of the nanoribbon; the rectangle marks one unit cell (supercell) 
of the ribbon. 



131 ] . The electron- ion interaction was described by the projector augmented wave method, 



and the energy cutoff was set to 500 eV. The structures were fully optimized using the conju- 
gate gradient algorithm until the residual atomic forces to be smaller than 10 meV /A. The 
supercell method with periodic boundary conditions was adopted to model the graphene 
nanoribbons (GNR), with a vacuum layer larger than 15 A to guarantee that there was no 
interaction between the GNR images in the neighboring cells. The edge energy is calculated 
as E edge = (E ribbon — E atom )/2L, where E ribbon is the total energy of the graphene nanoribbon 
in the supercell, E atom is the energy per atom in a perfect graphene without edge, and L 
is the length of edge. The edge stress is calculated as a e d ge = o- xx /2L, where a xx is the 
diagonal component of supercell stress tensor in the x-direction (defined along the edge), 
which is calculated using the Nielsen-Martin algorithm [ijj]. All other components of stress 
tensor vanish. We also note that DFT is suitable for calculating ground-state properties of 
lattice energies and stresses, to which the non-local many-body effects are not expected to 
be important. 

Figure 1 shows the calculated edge energy and edge stress of GNR armchair edges as 
a function of ribbon width ranging from ~ 3.5 to ~ 48 A. One notices that both edge 
energy and edge stress oscillate with the increasing width having a period of 3 but out 
of phase with each other. The oscillations are originated from the quantum confinement 



effect, as also manifested in the similar oscillations of electron band structures 
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FIG. 2: The AFM and PM zigzag edge stresses and edge energies of graphene nanoribbons as a 
function of ribbon width. Inset: schematics of the nanoribbon; the rectangle marks one unit cell 
(supercell) of the ribbon. 

oscillation magnitude of edge energy decays quickly with the increasing width and converges 
to ~ 1.0 eV/A, which agrees well with the previous first-principles valuesjs]. In contrast, 
the oscillation magnitude of edge stress decays much slower with a mean value of ~ -1.45 
eV/A(using negative sign as convention for compressive stress). The much larger oscillation 
in edge stress than in edge energy is possibly caused by the fact that edge stress equals to 
the derivative of edge energy with respect to strain, so that stress is much more sensitive to 
the width-dependent quantum confinement effect. There is also a slight revival effect in the 
stress oscillations at ~ 40A width, whose origin is not clear and needs further study. 

Figure 2 shows the calculated edge energy and edge stress of GNR zigzag edg es as a 
function of ribbon width ranging from ~ 5.0 to ~ 85 A. In this case, both edge energy and 
edge stress show very weak width dependence and converges quickly, again consistent with 
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jl8|. However, the zigzag edge is known 
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their corresponding electronic-structure behavior 

to have an antiferromagnetic (AFM) ground state [l7j. The AFM edge energy is calcu 
to be ~ 1.2 eV/A, about 0.2 eV/A lower than the paramagnetic (PM) edge energy ay 
Thus, we have calculated the edge stress of both AFM and PM states for comparison. We 
found that spin polarization have a sizable effect reducing the compressive stress from ~ 
-0.7 eV/A in the PM edge to ~ -0.5 eV/A in the AFM edge. 

Our first-principles stress calculations confirm qualitatively the recent empirical-potential 
results |lo| that both edges are under compressive stress. However, there are also some 
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significant differences. Two quantum manifestations of edge stress stand out, which are 
absent from the empirical prediction. One is the quantum oscillation of armchair edge stress 
with the increasing nanoribbon width, and the other is the reduction of zigzag edge stress by 
spin polarization. The physical origin of edge energy and edge stress is associated with the 
formation of one dangling bond on each edge atom. The repulsive interaction between the 
dangling bonds is believed to be one origin for the 'compressive' edge stress. In addition, in 
the armchair edge, it is well-known [20] that the edge dimers form triple -C=C- bonds with 
a much shorter distance (~ 1.23 A according to our calculation) adding extra compressive 
stress to the edge; while in the zigzag edge, spin polarization further reduces the compressive 
stress. Therefore, quantitatively, the armchair edge has a much larger compressive stress 
(~ -1.45 eV/A) than the zigzag edge (~ -0.5 eV/A). In contrast, the empirical potentials 
predicted a smaller compressive stress in the armchair edge (~ -1.05 eV/A) than in the 



zigzag edge (~ -2.05 eV/A)(lo|. 



The quantum effects in edge stress will in turn modify the mechanical edge instability. 
The compressive edge stress means the edge has a tendency to stretch. If we apply a 
uniaxial in-plane strain to a nanoribbon along the edge direction, the strain energy can be 
calculated as [lo| 

E str = 2r e Le + E e Le 2 + \E s Ae 2 (1) 

Here, A is the ribbon area, L is the edge length, r e is the edge stress, E e is the ID 
edge elastic modulus in a 2D nanoribbon, in analogy to the 2D surface elastic modulus in 



a 3D nanonlm 2 If. and E s is the 2D sheet elastic modulus. Since r e is negative, for small 
enough tensional strain e (positive), the negative first term (linear to e) in Eq. (1) can 
always overcome the positive second and third terms (quadratic to e) to make E str negative. 
Consequently, the ribbon is unstable against a small amount of stretching along the edge 
direction. Fitting first-principles calculations, by manually deforming the sheet and ribbon 
along the edge direction, to equation (1), we obtained E s ~ 21.09 eV/A 2 , .E e (amchair)^ 3 
eV/A and ^(zigzag) ~ 24 eV/A with r e already calculated above directly (see Figs. 1 and 
2). Our E s value is in good agreement with the experiment [22j] and empirical simulation 
result [h|. But our E e values are notably different from the empirical results [lo|. 

Another effective way to stretch the edge of a 2D sheet is by out-of-plane edge twisting 
and warping motions, which are barrierless processes. For example, assuming a sinusoidal 
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FIG. 3: (a) Armchair edge ripple amplitude versus ribbon width for A = 5GA. Inset: Schematics 
of ripple formation along the armchair and zigzag edge for A = 50 A. (b) Armchair (red) and zigzag 
(blue) edge ripple amplitude as a function of A. Light blue band shows the typical range of thermal 
fluctuation. 



edge warping with displacement fi e = asin(27rx/A) of amplitude a and wavelength A, which 
decays exponentially into the sheet as e~ y l l (See inset of Fig. 3), where I is the decay length, 
Shenoy et al. have shown that minimization of strain energy leads to characteristic length 
scales of such warping instability as I ~ 0.23A and a p 



-\T e )/(1.37E b + U.8EJX). Using 
their empirical-potential values of r e , E e and E s , they estimated that the warping magnitude 
of armchair edg e is smaller than that of zigzag edge and both are larger than typical thermal 
fluctuations [lOj. 

Our first-principle predictions, however, are different in several ways. First, absent from 
empirical prediction, the quantum oscillation of r e of armchair edge gives rise to an oscillating 
armchair edge warping amplitude for given wavelength as a function of nanoribbon width, 
as shown in Fig. 3a. Second, opposite to empirical prediction, the warping amplitude of 



(a) 



> 2 



Si 




edge stress 
edge energy 



1.2 



1.1 



1.0 



I 

13 



m 



20 40 60 80 

SW Concentration (%) 



0.9 



100 




FIG. 4: (a) The armchair edge stresses (with linear fit) and edge energies of graphene nanoribbons 
as a function of edge SW defect concentration, (b) The optimaized ribbon structure at the 50% 
SW defect concentration. 

armchair edge is much larger than that of zigzag edge, as shown in Figs. 3b. Third, the 
mechanical undulation of zigzag edges induced by compressive edge stress is comparable to 
thermal fluctuations (si, 0], as shown in Fig. 3b, and hence the two are difficult to distinguish. 
We also note that in addition to the continuum analysis we perform here, the accurate 
first-principles values of r e and E e can be used as input parameters for the finite element 
simulations of large graphene systems 10]. 

Because the compressive edge stress is partly originated from the dangling bond, natu- 
rally, we may saturate the dangling bonds to relieve the compressive stress. We have tested 
this idea by saturating the edge with H that indeed confirmed our physical intuition. For 
armchair edge in a 1-nm wide ribbon, we found H saturation changes the edge stress from 
-1.42 eV/A to -0.35 eV/A; for zigzag edge in a 2.0-nm wide ribbon, it changes the edge 
stress from -0.42 eV/A to +0.13 eV/A. Therefore, the H edge saturation, or saturation by 
other molecules in general, is expected to relieve the edge compression. Especially, it can 
reverse the compressive stress in a zigzag edge to tensile. 

Surface reconstruction has long been known as an effective mechanism in relieving surface 
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stress 



23j. Thus, in addition to H saturation, we have investigated possible edge reconstruc- 



tions in relieving the edge compressive stress. We have considered the Stone- Wales (SW) 



defect 



241 ]. which appealed to us become a SW defect in 2D is equivalent to a dislocation core 



in 3D that is known as a common stress relieve mechanism. Figure 4a shows the calculated 
armchair edge stress along with edge energy as a function of one type of SW defect (7-5-7 ring 
structure) concentration. Figure 4b shows an example of the optimized edge structure at the 
50% defect concentration. One sees from Fig. 4a that the edge stress increases linearly from 
compressive to tensile with the increasing SW defect concentration. The most stable edge 
structure is at ~ 25% defect concentration where the edge stress is very small and slightly 
compressive. A small stress value indicates that this chemically stable edge structure (with 
the lowest edge energy) is also most mechanically stable against deformation. 

Figure 5a shows the ground-state AF zigzag edge stress along with edge energy as a 
function of another type of SW defect (5-7 ring structure) concentration. Figure 5b shows 
an example of the optimized edge structure and spin charge density at the 50% defect 
concentration. Again, the edge stress increases linearly from compressive to tensile with the 
increasing defect concentration, the same as the case of armchair edge (Fig. 4a), but the 
edge energy decreases monotonically with the most stable edge having 100% of defects, in 
agreement with a recent first-principles calculation jg]. The initial compressive edge stress 
(~ -0.5 eV/A) is completely reversed to a large tensile value of ~1.2 eV/A in the most 
stable edge. Also, the 100% defected edge becomes non-spin-polarized. In general, the zigzag 
edge spin decreases continuously with the increasing SW defect concentration, similar to the 
behavior found previously for other types of defects 25]. 

We note that there is a first-principles calculation of graphene edge stress was reported 
recently 26] , but that calculation appears inconsistent with all other existing first-principles 
calculations . Their edge energies differ from all previous results [§, 9] and ours that are 
consistent with each other. The well-known spin-polarization of the zigzag edge was not 
considered, the quantum oscillation of armchair edge energy and stress was not shown, and 
their edge stress were not calculated from direct quantum-mechanical formula jl4l ]. 

In conclusion, quantum effects have been widely shown for electronic structure and ener- 
getic quantities of low-dimensional nanostructures. We demonstrate, in addition, quantum 
manifestations of mechanical quantities in graphene edge stress. We show that quantum 
confinement can lead to stress oscillations and spin polarization can reduce stress, which in 
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FIG. 5: (a) The zigzag edge stresses (with linear fit) and edge energies of graphene nanoribbons as 
a function of SW defect concentration, (b) The optimized ribbon structure and spatial distribution 
of spin density (charge density difference between spin-up and spin-down states in units of [13 A -2 ) 
of the AFM ground state at the 50% SW defect concentration. 

turn "quantum mechanically" modify the edge twisting and warping instability. We further 
show that H edge saturation and SW edge reconstruction can not only improve the 'chemi- 



cal' stability of graphene edges by lowering the edge energy as shown before|8j, |9[, but also 
enhance their 'mechanical' stability by converting compressive edge stress towards tensile 
and hence stabilizing the planar edge structure. Our first-principles findings, which cannot 
be captured by classical methods, provide new insights to the understanding of mechanical 
stability of graphene. We suggest that experimental measurement of armchair edge ripple 
magnitudes for different widths and comparison with zigzag edge ripples could confirm our 
predictions of quantum effects in graphene edge stresses differentiating from empirical re- 
sults. We expect the quantum manifestation of mechanical properties such as stress to exist 
generally in many low- dimensional nanostructures. 
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